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Recently, Holz and Wald |l| have presented a method 
for determining gravitational lensing effects in inhomogeneous 
universes. Their use of realistic galaxy models has been lim- 
ited to the singular, truncated isothermal sphere with a fixed 
mass. In this paper, their method is generalized to allow 
for matter distributions more accurately describing the ac- 
tual properties of galaxies, as derived from observations and 
N-body simulations. This includes the density profile pro- 
posed by Navarro, Frenk and White, as well as a distribution 
of galaxy masses. As an example of the possible applications 
of the method, we consider lensing effects on supernova lumi- 
nosity distributions. 

We find that results for different mass distributions of 
smooth dark matter halos are very similar, making lensing 
effects predictable for a broad range of halo profiles. We also 
note, in agreement with other investigations, that one should 
be able to discriminate smooth halos from a dominant com- 
ponent of dark matter in compact objects. For instance, a 
sample of 100 supernovae at redshift z — 1 can, with 99 % 
certainty, discriminate the case where all matter is in compact 
objects from the case where matter is in smooth halos. 
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I. INTRODUCTION 

Gravitational lensing has become an increasingly im- 
portant tool in astrophysics and cosmology. In partic- 
ular, the effects of lensing has to be taken into account 
when studying sources at high redshifts. In an inhomoge- 
neous universe, sources may be magnified or demagnified 
with respect to the case of a homogeneous universe with 
the same average energy density. 

The effects of gravitational lensing has been studied 
numerically by a number of authors, see, e.g., [0-0] ■ The 
most common method traces light rays through inhomo- 
geneous matter distributions obtained from N-body sim- 
ulations. Lensing effects are accounted for by projecting 
matter onto lens planes, and using the thin- lens approx- 
imation (see, e.g., |^]). 

Recently, Holz and Wald (HW; have proposed an- 
other ray-tracing method for examining lensing effects 



* E-mail address: goliath@physto.se 
-mail address: edvard@physto.se 



in inhomogeneous universes. This method can be sum- 
marized as follows: First, a Friedmann-Lemaitre (FL) 
background geometry is selected. Inhomogeneities are 
accounted for by specifying matter distributions in cells 
with energy density equal to that of the underlying FL 
model. A light ray is traced backwards to the desired 
redshift by being sent through a series of cells, each time 
with a randomly selected impact parameter. After each 
cell, the FL background is used to update the scale fac- 
tor and expansion. By using Monte Carlo techniques to 
trace a large number of light rays, and by appropriate 
weighting statistics for the apparent luminosity of 
the source is obtained. 

The advantages with this method are that light rays 
are traced through a three-dimensional matter distribu- 
tion without projection onto lens planes, thus avoiding 
any assumptions regarding the accuracy of the thin-lens 
approximation Furthermore, the method is flexible 
in the sense that cells may be taken to represent both 
galaxies and larger structures with different matter distri- 
butions, including non-spherical ones. For instance, HW 
have performed a number of tests to determine effects 
of clustering, and argue that this does not significantly 
affect statistical properties of magnification. They also 
investigate the case of substructure in the form of com- 
pact objects, and conclude that this can be adequately 
modelled by randomly distributed compact objects of ar- 
bitrary mass. It should be pointed out that the method 
is not well-suited to model clustering on scales larger 
than cell sizes. Still, galaxy clusters can be modelled by 
specifying appropriate masses with corresponding larger 
cells. Another drawback is that the method only con- 
siders infinitesimal ray bundles, making it impossible to 
keep track of multiple images. However, it is still possible 
to distinguish between primary images and images that 
have gone through one or several caustics jl]J| . 

HW considered pressure-less models with a cosmolog- 
ical constant, using the following matter distributions: 
point masses; singular, truncated isothermal spheres 
(SIS); uniform spheres; and uniform cylinders. The indi- 
vidual masses were determined from the underlying FL 
model using a fixed co-moving cell radius of R c = 2 Mpc, 
reflecting typical galaxy-galaxy separation length-scales. 

The aim of this paper is to allow for matter distribu- 
tions more accurately describing the actual properties of 
galaxies. We will extend the list of matter distributions 
to include the density profile proposed by Navarro, Frenk 
and White (NFW; §]) and we will use a distribution of 
galaxy masses. Also, other matter distribution parame- 
ters such as the scale radius of the NFW halo and the 
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cut-off radius of the SIS halo will be determined from dis- 
tributions reflecting real galaxy properties. The method 
of HW has also been generalized in Bergstrom et al. || to 
allow for general perfect fluids with non-vanishing pres- 
sure. 

Gravitational lensing effects may be of importance 
when, e.g., trying to determine cosmological parame- 
ters using observations of supernovae at high redshifts 
P, [lO|]I]| ] . In this paper, we study the effect from lensing 
on the luminosity distribution of a large sample of Type 
la supernovae at redshift z = 1. 



II. MASS DISTRIBUTION 

Realistic modelling of galaxies calls for realistic mass 
distributions and number densities, i.e., one has to allow 
for the possibility of the cell radius, R c , to reflect the 
actual distances between galaxies. 

An advantage of the method of HW is that it is very 
easy to allow for any mass distribution and number den- 
sity, including possible redshift dependencies, as long as 
the average density agrees with the underlying FL- model. 

Thus, for each cell we obtain a random mass, M, from a 
galaxy mass distribution dn/dM, and calculate the cor- 
responding radius from the condition that the average 
energy density in the cell should be equal to the average 
matter density of the universe at the redshift of the cell: 



47T , 
M = yfWcritiic, 



(1) 



where VLm is the normalized matter density, and p cr \t — 
3H 2 /8w is the critical density. A galaxy mass distri- 
bution can be obtained, for example, by combining the 
Schechter luminosity function (see, e.g., Peebles (L2J, 
Eq. 5.129) 



dn = 4>*y a e y dy, 
L 



(2) 
(3) 



with the mass-to-luminosity ratio (see, e.g., Peebles [ p.2| , 
Eq. 3.39) normalized to a "characteristic" galaxy with 
L = and M = M*. 
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Using Eq. (g), we find that 
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Assuming that the entire mass of the universe resides in 
galaxy halos we can write 



n(y)M(y)dy = p„ 



(7) 



2/min 



Using the Schechter luminosity function and the mass- 
to-luminosity fraction we get 



QmPc 



(8) 



Thus, by supplying values for n.*, reasonably well- 
determined by observations, y mm and y max , from which 
the dependence of M* is weak, together with parameters 
a and (3 we can obtain a Af* consistent with Qm- For 
the parameter values used in this paper (see Sec. [v|), we 
get 



M* « 7.50m • 1O 13 M . 



(9) 



III. THE NAVARRO-FRENK- WHITE 
DISTRIBUTION 

In the work of HW, the treatment of realistic galaxy 
models has been limited to the use of the singular, trun- 
cated isothermal sphere (SIS). Another often- used matter 
distribution is the one based on the results of detailed 
N-body simulations of structure formation by Navarro, 
Frenk and White @. The NFW density profile is given 
by 



p(r) 



Pcrit^e 



(r/R s )[l + (r/R s )Y 



(10) 



where S c is a dimensionless density parameter and R s 
is a characteristic radius. The potential for this density 
profile is given by 



$(r) = -ATT Pcrit 5 c R] 



,ln(l + x) 



const. 



(11) 



where x — r/R s . The matrix J a p, describing the evolu- 
tion of a light beam passing through a cell [see Eq. (37) 
in HW], can then be obtained analytically, see The 
mass inside radius r of a NFW halo is given by 



M(r) = 4Trp crit 6 c R 3 s 



ln(l 



1 



(12) 



Combining this expression with Eq. (jl]), i.e., setting M — 
M(r), we obtain 



9. 



M 



ln(l + x c ) 



l+x c 



(13) 



where x c = R c /R s . That is, for a given mass M , 5 C is a 
function of R s . From the numerical simulations of NFW 
we also get a relation between 8 C and R s . This relation 
is computed numerically by a slight modification of a 
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mag(filled beam) — mag(lensed) 

FIG. 1. Luminosity distributions for 10 000 perfect stan- 
dard candles at redshift z = 1 in a Q,m — 0.3, £Ia = 0.7 
universe. The magnification zero point is the luminosity in 
the corresponding homogeneous ( "filled-beam" ) model. The 
full line corresponds to the point-mass case; the dashed line 
is the distribution for SIS halos, and the dotted line is the 
NFW case. This plot can be compared with Fig. 22 of HW. 

Fortran routine kindly supplied by Julio Navarro. Of 
course, one wants to find a R a compatible with both the 
average density in each cell and the numerical simula- 
tions of NFW. Hence, we iteratively determine a value 
of R s consistent with both expressions for S c . Generally, 
R s will be a function of mass M, the Hubble parame- 
ter h, the density parameters £Im, ^a, an d the redshift 
z. However, we will use the result from Del Popolo p"3j ] 
and Bullock et al. pj that R s is approximately constant 
with redshift. We will compute R s for a variety of M, h 
and VLm (all at z = 0) in both open and flat cosmologies 
and interpolate between these values to obtain R s for any 
combination of parameter values. 



IV. TRUNCATION RADII FOR SIS-LENSES 

In their calculations for SIS halos, HW use a fix trun- 
cation radius d. However, using a realistic mass distribu- 
tion, the cut-off should depend on the mass of the galaxy. 
Here we derive an expression for d. 

The SIS density profile is given by 



psis(r) 



fill 

2tt r 2 ' 



(14) 



where a is the line-of-sight velocity dispersion of the mass 
particles. The mass of a SIS halo truncated at radius d 
is then given by 
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mag(filled beam) — mag(lensed) 

FIG. 2. Luminosity distributions for 10 000 sources at red- 
shift 2 = 1 in a Qm = 0.3, fl\ = 0.7 universe. This is the same 
situation as depicted in Fig. |l[ only that we have added an 
intrinsic luminosity dispersion of the sources with er m = 0.16 
mag. (corresponding to the case of Type la supernovae). 



M(d) = [ p(r)dV = 2a 2 d. 
Jo 



(15) 



We want this to be equal to the mass M given by the 
Schechter distribution, 



la d = M 



2ai \mJ \ a, 



(16) 



where we, in addition to M*, have introduced a char- 
acteristic velocity dispersion <r*. Combining the Faber- 
Jackson relation 



— = V 



with the mass-to-luminosity ratio, Eq. 
stitute for a in Eq. (|l6|), and obtain 



(17) 



we can sub- 



2ol 



M 



l-2A(l-/3) 



(18) 



Using Eq. (^), we can write the truncation radius for a 
halo with mass M = as 



d w 3.3fi A /Mpc. 



V. RESULTS 
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As an application of the method, we investigate lensing 
effects on observations of distant supernovae. In Fig. [I], 
we compare the luminosity distributions obtained with 
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FIG. 3. Luminosity distributions for 10 000 perfect stan- 
dard candles at redshift z — 1 in a Q,m = 1, ^a = universe. 
The magnification zero point is the luminosity in the corre- 
sponding homogeneous ( "filled-beam" ) model. The full line 
corresponds to the point-mass case; the dashed line is the dis- 
tribution for SIS halos, and the dotted line is the NFW case. 
This plot can be compared with Figs. 18 and 20 of HW. 

point masses, SIS lenses and NFW matter distributions 
in a Qm = 0.3, f^A = 0.7 universe, currently favoured 
by Type la supernova measurements [|l5[|l6| l . Sources are 
assumed to be perfect standard candles. The magnifica- 
tion, given in magnitudes, has its zero point at the filled 
beam value, i.e., the value one would get in a homoge- 
neous universe. Note that this value is cosmology de- 
pendent. Note also that negative values corresponds to 
demagnifications and positive values to magnifications. 
The point mass case (full line) agrees well with the re- 
sults of Holz and Wald. The SIS case (dashed line) is 
shifted towards the filled beam value when comparing to 
HW. This is due to the fact that the cut-off radii com- 
puted according to Eq. dig ) generally is much larger than 
the fix value of d = 200 kpc used by HW. Note that re- 
sults for SIS halos and NFW halos are very similar, even 
when we have no intrinsic luminosity dispersion of the 
sources. 

In Fig. ^| we have added an intrinsic luminosity disper- 
sion represented by a Gaussian distribution with a m — 
0.16 mag., due to the fact that Type la supernovae are 
not perfect standard candles. The effect is to make the 
characteristics of the luminosity distributions even less 
pronounced, since the form of the resulting luminosity 
distributions predominantly is determined by the form 
of the intrinsic luminosity distribution. It is still possible 
to observationally distinguish whether lenses consist of 
compact objects or smooth galaxy halos, as has been 
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FIG. 4. Luminosity distributions for 10 000 sources at red- 
shift z = 1 in a Qm = 1, ^a = universe. This is the same 
situation as depicted in Fig. |ij only that we have added an 
intrinsic luminosity dispersion of the sources with ti m — 0.16 
mag. (corresponding to the case of Type la supernovae). 



pointed out in |L7 18 1. Generating several samples con- 
taining 100 supernova events at z = 1 in a Qm = 0.3, 
SIa = 0.7 cosmology filled with smooth galaxy halos, we 
find that for 98 % of the samples one can rule out a 
point-mass distribution with a 99 % confidence level^. 
Furthermore, for a similar sample containing 200 super- 
novae, the confidence level is increased to 99.99 %. 

We have performed simulations for various cosmolo- 
gies, and found a substantial difference between SIS ha- 
los and NFW halos only in a matter-dominated universe, 
flM = 1, £7a = 0, where the luminosity distribution us- 
ing NFW halos is shifted towards the point-mass case 
(see Fig. |||). However, adding an intrinsic source disper- 
sion as in Fig. |J we see that even with the phenomenal 
statistics of 10 000 sources, it would be a difficult task 
to distinguish between the two density profiles. The in- 
creased number of high-magnification events for NFW 
halos would probably be the only way to make such a 
discrimination. 

In these calculations, we have used the following pa- 
rameter values (see further ||): 

• P = 0.2 



• a 



-0.7 



1 However, in 1 % of the samples, we will erroneously rule 
out the halo distribution with the same confidence level. 
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• 2/min = 0.5 

• 2/max = 2.0 



• n* = 1.9 • W~ 2 h 3 Mpc~ 3 

• <7* = 220 km/s 



• A = 0.25 

A more extensive discussion of the luminosity distri- 
butions of perfect standard candles obtained with the 
different halo models at different source redshifts can be 
found in Bergstrom et al. pi, where also some analyti- 
cal fitting formulas for the probability distributions are 
given. 



VI. DISCUSSION 

In this paper, the method of Holz and Wald [EJ has 
been generalized to allow for matter distributions reflect- 
ing the actual properties of galaxies, including the density 
profile proposed by Navarro, Frenk and White Q . In or- 
der to make matter distributions as realistic as possible, 
all parameter values in the lens models are obtained from 
reasonable probability distributions, as derived from ob- 
servations and N-body simulations. This includes the 
mass of the galaxies, the truncation radius of SIS lenses 
and the characteristic radius of NFW halos. One of the 
virtues of this method is that it can be continuously re- 
fined as one gains more information about the matter 
distribution in the universe from observations. 

The motivation for these generalizations is to use this 
method as part of a model for simulation of high-redshift 
supernova observations. In this paper, we have con- 
sidered lensing effects on supernova luminosity distribu- 
tions. Results for different mass distributions in smooth 
dark matter halos was found to be very similar, making 
lensing effects predictable for a broad range of density 
profiles. Furthermore, given a sample of 100 supernovae 
at z ~ 1, one should be able to discriminate between the 
case with smooth dark matter halos and the (unlikely) 
case of having a dominant component of dark matter in 
point-like objects. 
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